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ABSTRACT 


We present distributions of two molecular clouds having velocities of 2kms“^ 
and 14kms“^ toward RCW 38, the youngest super star cluster in the Milky 
Way, in the J =1-0 and 3-2 and J =1-0 transitions. The two 

clouds are likely physically associated with the cluster as verihed by the high 
intensity ratio of the J=3-2 emission to the J=l-0 emission, the bridging feature 
connecting the two clouds in velocity and their morphological correspondence 
with the infrared dust emission. The velocity difference is too large for the 
clouds to be gravitationally bound. We frame a hypothesis that the two clouds 
are colliding with each other by chance to trigger formation of the ~20 O stars 
which are localized within ~0.5pc of the cluster center in the 2kms“^ cloud. 
We suggest that the collision is currently continuing toward part of the 2kms“^ 
cloud where the bridging feature is localized. This is the third super star cluster 
alongside of Westerlund2 and NGC3603 where cloud-cloud collision triggered 
the cluster formation. RCW 38 is the youngest super star cluster in the Milky 
way, holding a possible sign of on-going O star formation, and is a promising site 
where we may be able to witness the moment of 0-star formation. 

Subject headings: ISM: clouds — radiolines: ISM 
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Introduction 


1.1. Important role of high-mass stars 


High-mass stars are very energetic and influential on dynamics of the interstellar 
medium via stellar winds, ultraviolet radiation and supernova explosions. Injection of the 
heavy elements at the end of their lives is another important effect of high-mass stars 
in galactic evolution. It is therefore one of the most important issues in astrophysics to 
understand formation of high-mass stars and considera 
numb er of works published so far (for recent reviews see 
201411 . 


3le efforts have been mac 

: e bv a 

Zinnecker & Yorke 

2007; 

Tan et al. 


1.2. Theoretical attempts for understanding high-mass star formation 


It is now generally accepted that most of the high-mass stars are formed by mass 


accretion instead of collisional s 


paradigm, monolith ic collapse flYorke fc Sonnhalter 


flBonnell et al. 


ellar merging after some 


years of debates. In the current 


20021) and competitive accretion 


20041) are the two scenarios of the mass accretion processes. In the 


monolithic collapse a massive dense molecular cloud is the initial condition for the star 
formation, which is followed by disk accretion, whereas in the competitive accretion a 
massive cluster of stars/pre-stellar cores is assumed to grow in mass by accretion of the 
ambient gas. The two scenarios still lack detailed confrontation with the initial conditions 
provided by observations. This is in part due to the lack of well-resolved observational 
views of high-mass star formation at small distances; even the Orion Nebular Cluster, the 
most outstanding high-mass star forming region at the smallest distance 400 pc, is three 
times more distant than typical low-mass star-forming regions like Taurus, and the second 
nearests are at a distance around 2kpc. Therefore, our understanding of high-mass star 
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formation yet remains to be significantly corroborated. 


It is interesting to note that insight-full disc ussion is made on t h e mec hanisms of 
high-mass star formation in the review article bv IZinnecker fc Yorkel (120071) . where we 


hnd the following sentence: “Rapid external shock compression (i.e., supersonic gas 
motions) generating high-column densities in less than a local free-fall time rather than 
slow quasi-static build-up of massive cores may be the recipe to set up the initial conditions 
for local and global bursts of massive star formation.” This idea of supersonic gas motions 
to collect gas was not explored further in the literature in the context of high-mass star 
formation except for the several papers on cloud-cloud collision discussed in Section 1.3. 


1.3. Increasing evidence for clond-cloud collision as a triggering of high-mass 

star formation 


Molecular observations with NANTEN2 have shown that 0-type star formation in 
the Milky Way disk is triggered by supersonic collisions between two clouds at velocity 
of 10-30 km s“^ in several places. The hrst discovery was made in the super star cluster 


West erlund 2 which is ionizing an Hll region ROW 49 flFnrnkawa et al. 


2009 


Ohama et al 


2010), and the authors suggested that two giant molecular clouds collided to trigger 
formation of the cluster some 2 Myrs ago. These observations were followed by discovery 
of cloud-cloud collisions which t riggered formation of the super star cluster/mini-starburst 


in NGC3603 (Fukui et al. 


(iTorii et al. 


3, 


2014 , and forma 


2011 ) and RCW120 ( Torii et al. 


ion o f the single O star in Hll regions M20 


20151) . Most recently, in the Large Magellanic 


Cloud ALMA observations led to a discovery of formation of a 40 Mq star N159W-S 
triggered by collis ion between two hlamentary clouds at a relative velocity of ~ 10kms“^ 


flFukui et al 


20151) . These results raise a novel possibility that formation of O stars may be 


due to, or at least facilitated, by the supersonic compression in cloud-cloud collision, which 
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was not considered in the theories for high-mass star formation, either in the monolithic 
collapse or in the competitive accretion scenarios. On a global scale of galaxies, nnmerical 


simnlations show tha 


(Tasker fc Tan 


2009; 


collisions betwee n clonds are fairly freqnent at every ~10Myrs 


Dobbs et al 


of gia nt molecular clouds flFnkni et ah 


2015), wh i ch is similar to an evo l utionary timescale 


1999 


Kawamnra et ah 


2009 


Fnkni fc Kawamnra 


20101 ). It is thus becoming possible to consider cloud-cloud collision as one of the modes of 
high-mass star formation. 


Cl oud-cloud collis i on wa s stud ied using hydrodynam ics simulations by 


(1199211 . 


Anathpindikal (1201011 . and 


Takahira et al. 


Habe&Ohta 


(1201411 . These authors found that 


cloud-cloud collision induces the formation of cloud cores by enhanced self-gravity as 
a consequence of shock compression, while the connection of cloud-cloud collision to 
high-mass star formation was not discussed into depth. Recent magnetohydrodynamical 
(MHD) numerical simulations on cloud-cloud collision show that colliding molecular gas 
can indeed create dense and massive cloud cores, precursors of high-mass stars, in the 
shock-compressed int erface, giving a theoretical basis for triggered the 0-star formation 


(jlnoue fc Fnkni 


20131) . In the simulations after the collision the density increases from 
300-1000 cm“^ to > 10^cm“^ and the turbulent velocity and magnetic held increase by 
about a factor of hve, leading to a 100-times higher mass-accretion rate to form stars. This 


indicates that the conditions for t 


formation theories (iMcKee fc Tan 


le hig h mass-accretion rate postulated in high-mass star 


20031) are created by cloud-cloud collision. 


Following these hndings, as a next step we need to have a better understanding on the 
physical details of the collision, which include direct observations of the shock interaction as 
well as the statistical properties of cloud-cloud collisions. We are carrying out a systematic 
study of high-mass star formation toward Hll regions in order to elucidate the role of 
cloud-cloud collision with NANTEN2, Mopra, and ASTE telescopes. These Hll regions 
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include Spitzer bubbles, compact Hll regions, and mini-starbursts including NGC6334, 
NGC6357, W43, etc. 


1.4. RCW 38 as the youngest super star cluster in the Milky Way 


Hoping to better understand the origin of massive clusters we in this paper focus 

le you ngest 


ourselves on the super star cluster RGW 38 (iRodgers et ah 


super star cluster with an age o: 


Table 2 in 


Portegies Zwart et al 


<1 M yr in the Milky Way ( WoIk et al 


196nh. whic 


n IS 


2006 


2008 


see also 


201011 . 


citetwol2006 estimated the age of RGW 38 as O.SMyrs with an upper limit of 1 Myrs so that 


i t doe s not conflict with the number of O stars for the size of the cluster, and 


Getman et al. 


(120141 1 measured the ages of the several members of the RGW 38 cluster, which ranges 


from 0.1-3.5 Myrs. The cluster age may be even smaller than 0.5 Myrs as discussed in the 
present paper (Section 4). The extreme youth will beneht study of the molecular gas close 
to the central O stars, which is particularly important to address the ionization/dispersal 
of the molecular gas. RGW 38 has large total mass, similar to that of the Orion Nebula 
Gluster, and is an active high-mass star forming region, second nearest to the sun, at a 
distance of 1.7kpc. An over view of the previous studies on RGW 38 is summarized in an 


article by 


Wolk et al 


(1200811 . 


RGW 38 is known as one o 


10,000 stars ( Kuhn et al. 


2015 


the closest 


Broos et al. 


ligh-mass star forming regions containing 


2013 


Wolk et al. 


2006 


; see also 


Lada fc Lada 


20031) (Figure!!]) . In the central part o f RGW 38, two remarkable infrared peaks have 


been identihed (iFrogel fc Persson 


1974 1. The brightest at 2/mi is labeled I RS 2, which 


corre s ponds to an 05.5 bin ary located at the center of the RGW 38 cluster (IDeRose et al. 


200911 . 


Furniss et al. 


(I 1975 I) derived a total infrared luminosity toward IRS 2 to be 7 x 10^ Lq. 


The brightest feature at 10 /mi is found 0.1 pc west of IRS 2, and is labelled IRS 1. IRS 1 is 





























































a dust ridge extending by 0.1-0.2pc in the north-south direction (Figure[I}3^ 
a col or temperature of about 175 K and includes several dust condensations 


, which has 


Smith et al, 


19991) . Infrared and millimeter wave observ ations has indicated a “ring-like” or “horseshoe 


struc ture around IRS 2 about 1' - 2' across (IHuchtmeier 


Huchtmeier 

1974; 

Vigil 

2004; 

Wolk et al. 

2006 


20081 ). Inside the ring-like shape two regions are cleared of dust and form cavities; one is the 


region centerec 


a similar size fjSmith et al 


on IRS 2 with a c 


1999 


iameter of ^0.1 p c, and another is just west of IRS 1 with 


Wolk et al 


20081) . suggesting that they were formed by 


feedbacks from high-mass stars. A large-area inf rared image shows t hat RCW 38 includes 


numerous hlaments and bubbles (see Figured^). 


Kaneda et al. 


(120131) made comparisons of 


the dust and polycyclic aromatic hydrocarbon (PAH) maps with the C[II] 158/rm map and 
discussed that the PAHs and dust grains are well mixed along the line-of-sight and PAHs 
play an important role for photo-electric heating of gas in photo-dissociation regions. 

Many observational studies to investigate the cluster members of RCW 38 have been 
carried out so far. Near infrared observations by using the Very Large Teles cope (VLT) 


ident ihed more than 300 young stars for t he ~0.5pc^ a rea ce ntered on IRS 2 (IDeRose et al. 


2009 

). Chandra observations reported by 

Wolk et al. 

(2006) 


X-ray sources likely associated with RCW 38, and com pleteness arguments 


estimated cluster size of between 1500 and 2400 stars. 


Winston et al. 


suggest a total 


( 20111) identihed 624 


young stellar objects (YSOs) by the SpitzerlYRkC observations toward a 30' x 30' region 
of RCW 38. In a part of The Massive Young star-forming Complex Study in Infrared and 
X-rays (MYStlX), 886 pre-main sequence stars have been identified in RCW 38, wl 


concentratec 

. wit 

lin 

Kuhn et al. 

2013 

). 


identihed by 


lin 0.5 pc of the cluster center (IFeigelson et al. 


2013 


RCW 38, w 

rich £ 

Broos et al. 

2013 


Wolk et al. 


the high- mass stars, nearly 60 0-star candidates have been 


(120061) and 


Winston et al. 


(20111 ) in a large a rea of RCW 38. O ut 


of them, about 20 0-star candidates are found in the central ~0.5pc (jWolk et al 


20061 ). 


citetkuh2015 also pointed out that the RCW 38 cluster shows a particularly high central 















































































stellar density compared with the other nearby clusters with O stars (e.g., the Orion 
Nebula Cluster, W40, RCW36, M17, etc.). It seems natural that such an O star cluster 


has tens of O stars as cluster members, w 


star spectroscopically flDeRose et ah 


rile only IRS 2 is confirmed to be a double 05.5 


20091) . We note cautiously that the number of O 


stars is still to be corroborated observationally by considering the uncertainty of the stellar 
mass determinatio n based on pre-main-sequence models for the K-band stellar magnitude 


llWolk et al 


20061) 


Molecular line observations toward RCW 38 are very limited. 


Yamagnchi et al. 


fll999h 


observed a large area of RCW 38 in the ^^CO J=l-0 emission with NANTEN and found 
a molecular cloud at 2-3kms“^ which coincides with IRS2, although detailed structures 
including the ring-li ke component we r e not resolved due to a relatively coarse resolution. 


On the other hand. 


Zinchenko et al, 


(119951) observed the central region of RCW 38 at 
HPBW of ~ V in the CS emission to resolve the ring shape surrounding IRS 2. ^^CO 


J=l-0 observations with S 
by 


ST a t 45" resolution for a small 160" x 160" region performed 


Gynlbndaghian fc May! (120081 ) show two molecular clouds likely associated with RCW 38 
at two different velocity ranges of —3--|-2kms“^ and -|-3--|-8kms“^. 


In this paper, we present new molecular line observations toward RCW 38 obtained 
with NANTEN2, Mopra, and ASTE telescopes. The NANTEN2 CO J=l-0 observations 
covered a large area of RCW 38 at an angular resolution of 4' to describe the large 
scale molecular morphologies and velocities, while the Mopra CO J=l-0 and ASTE CO 
J=3-2 observations at angular resolutions of 22"-35" revealed the detailed molecular gas 
distributions and physical conditions such as temperature and density around IRS 1 and 
IRS 2. Section 2 describes observations and Section 3 observational results and the physical 
properties of the molecular clouds. Section 4 gives discussion on cloud-cloud collision and 
O star formation, and Section 5 concludes the paper. 
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2. Observations 

2.1. NANTEN2 CO J=l— 0 observations 

The NANTEN2 4-m mm/sub-mm telescope situated in Chile was used to observe 
a large area 1.8° x 1.8° of RCW 38 in the and J=l-0 transitions during a 

period from May 2012 to December 2012. The 4K cooled SIS mixer receiver equipped 
with NANTEN2 provided a typical system temperature Tgys of ~250K in DSB. The 
backend was a digital spectrometer which provided 16384 channels at 1 GHz bandwidth and 
61kHz resolution, corresponding to 2600 kms“^ and 0.17kms“^, respectively, at 110 GHz. 
The obtained data were smoothed to a velocity resolution of 0.8kms“^ and angular 
resolution of 240". The pointing accuracy was checked to be better than 15" with daily 
observations toward the Sun and IRC 10216 (R.A., Dec.) = (9*^47™57!406, 13°16'43"56). 
The equinox of the celestial coordinates used in this paper is J2000.0. The absolute intensity 
calibration was done with daily observations of the p-Ophiuchus molecular cloud (R.A., 
Dec.) = (16^32™23!3, —24°28'39"2) and Perseus molecular cloud (R.A., Dec.) = (3^29™ 19®0, 
31°24'49"0). 


2.2. Mopra CO J=l—0 observations 

Detailed CO J=l-0 distributions around RCW 38 were obtained by using the 22-m 
ATNF (Australia Telescope National Facility) Mopra mm telescope in Australia at angular 
resolution of 33" in July 2012. We simultaneously observed the ^^CO J=l-0, ^^CO J=l-0, 
and C^®0 J=l-0 transitions toward a 11' x 15' area of RCW 38 in the OTF mode with a 
unit held of 4' x 4'. Tgys was 400-600 K in the SSB. The Mopra backend system “MOPS” 
which provided 4096 channels across 137.5 MHz in each of the two orthogonal polarizations 
was used in the observations. The velocity resolution was 0.088 km s“^ and the velocity 


10 


coverage was 360kms“^ at 115 GHz. The obtained data were smoothed to a HPBW of 40" 
with a 2D Gaussian function and to 0.6kms“^ velocity resolution. The pointing accuracy 
was checked every 1 hour to keep within 7" with observations of 86 GHz SiO masers. We 
made daily observations of Orion-KL (R.A., Dec.) = (—5^35™14!5, —5°22'29'.'6) to estimate 
“extended beam efficiency” by making compariso ns with the peak temperature of 100 K in 


Orion-KL shown in Figure 6 of 
efficiency of 0.48. 


Ladd et ah 


(120051) . We hnally obtained an extended beam 


2.3. ASTE J=3—2 observations 


Observations of ROW 38 in the ^"^GO J=3-2 tr ansition were perfo r med 


the ASTE 10-m telescope located in Ghile in 2013 flEzawa et al 


2004 


by using 


20081) . The 


waveguide-type sideband-separating SIS mixer recei ver for the 


‘GAT345” having system temperature of ~ 250 K fjinoue et al. 


single 


sideband (SSB) 


20081) and the digital 


spectrometer “MAG” with the na.rrow- band mode providing 128 MHz bandwidth and 


20001 ) . which correspond to lllkms ^ velocity coverage 


0.125 MHz resolution (iSorai et al. 
and 0.11 km s“^ velocity resolution at 345 GHz, were used. The observations were made 
with the GTE mode at a grid spacing of 7.5", and the HPBW was 22" at the ^^GO J=3-2 
frequency. The pointing accuracy was checked every ~ 1.5 hours and kept within 5" with 
observations of RAFGL5254 (R.A., Dec.) = (09^13™53®94, —24°51'25'.'l). The absolute 
intensity calibration was done with observations of Orion-KL (R.A., Dec.) = (—5’^35™14!5, 
—5°22'29'.'6), and the day-to-day fluctuations of the peak intensity were within 10%. 
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3. Results 


3.1. Large-scale molecular distributions with NANTEN2 


We first present large-scale molecular distributions of RCW 38 using the NANTEN2 
J=l-0 dataset. Compared with the NANTEN observations given bv lYamaguchi et al. 


fll999l) . although the beam size of 180" is the same, a grid spacing of 60" in the present 


NANTEN2 observations is a third of the size, which enables us to reveal the gas distributions 
into more detail. 


As shown in the CO distributions in FigureEj toward the direction of RCW 38, two 
CO clouds are distributed at radial velocities ulsr of —4--|-8kms“^ and -|-9--|-14kms“h 
In Figure[2t^a), the blue-shifted cloud has a strong peak just coinciding with the central 
part of RCW 38. The compact peak with a size of 0?1 (~ 3 pc at RCW 38) has a feature 
elongated t oward the southwest of R CW 38, showing a “head-tail” structure, which was also 


detected in 


Yamaguchi et al. 


(119991 ). The head-tail structure is surrounded by the diffuse 
CO emission distributed at around ulsr of 0--l-8kms“\ covering the present observed 
region. 


Contrary to the blue-shifted cloud shown in Figure[2](a), the red-shifted cloud has weak 
and dispersed CO distribution. A CO peak at (/, b) ~ (268?0, —1?1) with relatively strong 
emission is located just at the southeast of RCW 38, pointing toward the center of RCW 38. 
This CO peak is well covered with the Mopra and ASTE observations as indicated by a box 
with solid lines in Figure[2])b). There are several other CO components around the CO peak, 
and they all seem to be continuously distributed along the east-west direction roughly at a 

range of {l,b) ~ (267?5 -268?4, —1?0). In the Dec.-velocity diagram in Figure[2](c), the 

red-shifted cloud can be identihed separately from, although it looks being connected with, 
the blue-shifted cloud and the diffuse emission, suggesting a possible physical relationship 
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of the red-shifted cloud with RCW 38. 


3.2. Detailed molecular distributions with Mopra and ASTE 


We present the detailed molecular distributions around RCW 38 using the Mopra and 
ASTE datasets. In Figure[3l the integrated intensity distributions (left and central panels) 
and the Dec.-velocity diagram (right panels) are shown for the J=l-0 emission 

(upper panels) and the J=l-0 emission (lower panels). The R.A.-velocity diagram 

in the J=l-0 emission is also shown in FigurelH In Figures[3](a) and[3](d), the 

outstanding compact peak seen in Figure[2](a) is resolved. It shows a ring-like structure 
having a diameter of 1-2pc, with a cavity centered on IRS2, and several hlamentary and 
bubble-like structures are radially extended outside the ring-like structure by l-2pc. In the 
position-velocity diagrams in Figures[3](c) and HI intensity decrease is seen at around 
ulsr of 2-4kms“^. The emission in Figure[31(f) shows a single peak at ~ 2kms“^, 
indicating that the dip at the same velocity in the is due to self-absorptioiji] The 

blue-shifted cloud has a peak at 2kms“^. The ring-like structure in consists of four 

clumps (clumps I-IV), and the ring-like structure itself has an elongated shape along the 
northeast-southwest direction (Figure H]). The size of the clumps ranges from 0.3 pc to 
0.5 pc, and the linewidth An from 4kms“^ to 5kms“^ (Table [T]). 


In contrast to the complicated distributions of the 2kms“^ cloud, the 14kms“^ cloud 
has a simple ridge elongated from the north to the south (Figures|3]D and|3^). We shall 
hereafter call the 2kms“^ cloud and the 14kms“^ cloud as the ring cloud and the finger 
cloud, respectively, by their shapes. The hnger cloud is not likely affected by self-absorption 


^The two velocity features reported by iGvulbudaghian fc Mavl (120081) in J=l-0 


correspond to the two peaks caused by self-absorption. 
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in ^^CO, since and show similar distribntions in the Dec.-velocity maps in 

FignresEI^c) andEI^f). We hnd three “bridging” featnres in velocity between the ring cloud 
and the huger cloud at Dec. ~—47° 30f0 (toward the cluster), —47° 34^5 and —47° 3915, 
respectively, as indicated by arrows in Figure [3](c) and they are not seen beyond the velocity 
of the huger cloud. Their positions are denoted as BRl, BR2 and BR3 in Figure El The 
bridging feature distributed toward the cluster is also seen in the R.A.-velocity diagram in 
FigurelH At its bottom, the ridge intersects another elongated structure running from the 
east to the west, and these two orthogonal elongated structures comprise the huger cloud 
in the present observed area (FiguresEl(b) and (e)). As shown in FiguresEl(c) andE](f), the 
huger cloud shows a uniform velocity gradient of ~1 kms“^ pc“^ along the ridge. 


3.3. Intensity ratios of the J=3—2 emission to the J=l—0 emission 


In order to investigate the physical properties of the ring cloud and the huger cloud, we 
present distributions of the intensity ratio of the J=3-2 emission to the J=l-0 


useful diagnose o f the molecular gas properties (e.g.. 


Fukui et al. 


Oha,ma et al. 

201 Q; 

Torii et al. 

2011 


201411 . 


As seen in the Dec.-velocity diagram of i? 3 _ 2 /i-o in Figure[7l both of the ring cloud 
and the huger cloud have typical ratios of 0.6-0.8, up to over 1.0. Since the ring cloud is 
strongly ahected by the self-absorption in spectra around its central velocity range, 

~0-4kms“^, the ratio in this velocity range is not reliable. On the other hand, low 
intensity ratios are seen in the northern part (Dec. > —47° 2710) and the southern part 
(Dec. < —47°34l0) of the ring cloud, where the dihuse CO emission is widely distributed, 
while the huger cloud retains its relatively high intensity ratios of ~0.6-0.8 throughout the 
cloud. 
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3.4. Temperature and density of the molecular gas 


We utilize the large velocity gradient (LVG) analysis (e.g., 


Goldreich fc Kwan 


1974h 


to estimate kinetic temperature and molecular number density n(H 2 ) of the two clouds 
in order to interpret their intensity ratio distributions. We hrst present curves of i? 3 _ 2 /i_o 
as a function of Tk in Figure[ 8 ]for various densities from lO^crn”^ to lO'^cm”^, where 
X(CO)/{dv/dr) = 10“^ fkm s"^ pc~^)~^ is assurned. We adopt an a bundance ratio X(GO) 


= |'^CO|/|H 2 l = 10-1 (e.g., 


Frerking et ah 


1982 


Leung et ah 


19841) . and typical dv/dr is 


estimated for clumps I-IV as 4kms ^/0.4pc = 10kms ^pc ^ 


FigureE] provides a guide for the i? 3 _ 2 /i-o distributions in Figures[71 If i? 3 _ 2 /i-o is 
larger than 0.7, which is depicted by a dashed line in FigureEl 7k is always higher than 
10 K for every n(H 2 ). 10 K is a canonical value of the galactic molecular gas without 
signihcant star formation, suggesting that the gas with Tk higher than 10 K is caused by 
some additional heating. On the other hand, if R 3 - 2 / 1-0 is as small as 0.4, Tk cannot be 
determined, since it has a large variation from <10 K to >100 K depending on n(H 2 ). Both 
of the ring cloud and the huger cloud show R 3 - 2 / 1-0 higher than 0.7 up to over 1.0 at many 
points, suggesting that these two clouds have Tk higher than 10 K. It is reasonable that 
the heating is mainly due to the O stars in RGW 38. On the other hand, the diffuse GO 
emission surrounding the ring cloud with i? 3 _ 2 /i_o of 0.4 is of low excitation. 


In order to provide more quantitative constraints for Tk and ? 7 ,(H 2 ), we add the ratio 
of ^^GO J=l-0 to ^^GO J=l-0 (hereafter T 1 3 / 19 ) in the LVG analy sis, with an assumption 


of the abundance ratio [^^GO]/[^^GO] = 77 ( Wilson fc RoodI 11994 ). Ten target regions 


A-J for Tk and n(H 2 ) estimates in the two clouds. The GO spectra and the LVG results 
are presented in Figuresl9] and ITO] for the individual target regions. In the GO spectra, the 
velocity ranges used for the LVG analysis are shown by shade. For the ring cloud, the hve 
target regions are chosen to have apparent self-absorption in their ^^GO spectra. For the 
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ring cloud, regions A-C are taken at the high i? 3 _ 2 /i_o regions around the rim of the dense 
part, while regions D and E are at the low -R 3 _ 2 /i-o regions widely distributed in the north 
and the south. Regions F-J in FigurefTUlare taken to cover a large area of the huger cloud. 
In the diagrams of the LVG results, R 3 _ 2 /i-o and R 13/12 distributions are shown with the 
blue lines and red lines, respectively, and their errors are shown with the colored area. 
The errors are estimated with 1 a noise huctuations of the spectra and the 10 % relative 
calibration error for each CO transition. Since the J=l-0 emission and the 

J=l-0 emission were taken simultaneously with the same receiver and the same backend at 
Mopra, we assume that the 10% relative calibration error is canceled for R 13/12 and adopt 
the error only for Rs- 2 /i-o- In addition, dv/dr is estimated for the individual regions, and 
X(CO) is assumed to be 10“^ same as in Figured Finally, Tk and n(H 2 ) are given at the 
region where R 3 _ 2 /i-o curve and R 13/12 curve overlap. 

The results are summarized as follows: Regions A-C in the ring cloud and all the 
regions in the huger cloud (regions F-J) show Tk higher than 10 K, typically 30-40K, up 
to more than 50 K. Typical n(H 2 ) for these regions are about 10^-10^cm“^, where only 
lower limits are given for regions A, B, and F2 and no solutions are given for regions FI 
and I. Regions D and E located in the south of the ring cloud show signihcantly low Tk less 
than 10 K and n(H 2 ) of ~ 3 x 10^cm“^. 

As a summary, both of the ring cloud and the huger cloud show remarkably high 
temperatures of Tk > 20-30K. Some heating source/mechanism is necessary to understand 
the results, and radiative heating by the O stars in RCW 38 is a reasonable mechanisum. 
The present results suggest that, despite of a large velocity separation of ~12kms“^, both 
of the ring cloud and the huger cloud are associated with RCW 38. 
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3.5. Molecular masses of the ring cloud and the finger cloud 


By adopting the distance of RCW 38, 1.7kpc, to both the ring clond and the finger 
clond, their molecnlar masses are estimated. For a large scale, the NANTEN2 J=l-0 
data presented in Figures[2](a) and [2Jb) are used to estimate the total mole cular masses 


of th e two clouds. We use an X(CO) factor of 2 x 10^°cm ^ (Kkms ^ ( Strong et ah 


19881 ). which is an empirical conversion factor from J=l-0 integrated intensity into 


H 2 column density. The masses estimated with the X(CO) factor (hereafter Mx(co)) of fhe 
head-tail structure of the ring cloud and the finger cloud are calculated as 3.0 x 10^ Mq 
and 1.9 X 10^ Mq toward the regions enclosed with dashed lines in Figure[2](a) and[2](b), 
respectively. Mx(co) of the finger cloud is more than one order of magnitude smaller than 
that of the ring cloud. 


The molecular masses of the central region of RCW 38 shown in Figure[3](b) are 
estimated with the Mopra CO data set. Mx{co) of the finger cloud is estimated to be 
1.2 X IO^Mq at ^^CO J=l-0 integrated intensity hF(^^CO 1-0) > 10Kkms“^ (= 3cr). 
On the other hand, because the ^^CO J=l-0 intensity is reduced by by the strong 
self-absorption particularly in the central l-2pc of the ring cloud (see Figure[31(c)), we 
estimate the mass of the ring cloud using the Mopra ^^CO J=l-0 data assuming the 
local thermodynamic equilibrium (LTE). In estimating the molecular mass with the LTE 
assumption (hereafter Mlte), excitation temperature Tex must be provided, and the peak 
intensity of the ^^CO J=l-0 emission T(^^CO) is usually used to estimate Tgx with the 
following equation. 


ln{l + 5.53/(T(i2CO) +0.819)}' 

Despite of the strong self-absorption in the ^^CO profiles, observed brightness temperature 
in the ^^CO clumps is typically very high up to 45 K in clump III, which corresponds to Tex 
of 49 K in equation ([1]). This figure is roughly consistent with the LVG analysis, resulting 
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in Tk > 30-50K (Figure[9]) in regions A-C which are located in the ontskirts of the ring 
cloud. We therefore adopt a uniform Tex of 49 K and estimate Mlte of the dense ring 
structure of the ring cloud in Figure[3](d) to be 5.2 x 10^ Mq at J=l-0 integrated 

intensity hF(^^COl-O) >22.5Kkms“^. Mlte of the four clumps in Figure[5]^a) are 
estimated to be 1.0 x 10^ Mq for clump I, 8.8 x 10^ Mq for clump II, 5.7 x 10^ Mq for clump 
III, and 5.6 x 10 ^Mq for clump IV at hF(^^COl-O) >32.5Kkms“^ (see also Tabled]). The 
H 2 column density V(H 2 ) of the clumps is typically as high as 1 x 10^^cm“^. 


3.6. Comparisons with infrared images 

Figured!] shows a comparison between the J=l-0 distribution and the Spitzer 

3.6 /im image which consists of several hlamentary features around the central peak. Only 
panel (c) is for ^^CO. Figure fTTf a). a large-scale view, shows that the ring cloud and the 
intense 3.6 /rm emission corresponds well with each other. The huger cloud in FigredUb) 
shows possible correspondence with the southern part of the 3.6/im distribution at (R.A., 
Dec.) = (08^ 59™ 5®, —47° 35'). Figure fTlT cl compares the YSOs and O star candidates in the 
ring cloud. The YSO distribution corresponds to the overall distribution of the ring cloud, 
and the O stars with the northwestern part of the cavity of the ring cloud. Figure dU dl 
indicates that nearly 20 of the O stars are located toward the tip of the huger cloud 
overlapping with the bridging feature, while the distribution of the O stars shows an ohset 
slightly toward the south of the peak of the huger cloud. The other YSOs do not show 
strong correlation with the huger cloud in the south iFigure fTlT cT. 

FignrefT^ shows velocity channel distributions of CO overlaid on the infrared image. 
We label the nine infrared hlamentary features by a-h in the hrst panel of Figured^] The 
CO features generally show good correspondence with the infrared features as listed in 
Table 2. Correspondence of the infrared image with the huger cloud is not so obvious 
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in Figure[T2l FigurefT^ an overlay of the ring cloud and the finger cloud with the VLT 
image at a smaller scale, however, shows a good agreement between the finger cloud and 
IRS 1 which is elongated in the north to the south while the angular resolution of the 
CO distribution is too coarse for spatially resolving the bridging feature to a resolution 
comparable to the infrared image. We interpret that the dark lane toward IRS 1 along with 
another dust feature comprising the western cavity corresponds to the finger cloud which 
probably lies on the near side of the cluster. IRS 1 is clearly irradiated by the cluster and 
the physical association of the bridging feature with the cluster is supported. 


Discussion 


4.1. RCW 38 as the youngest snper star cluster in the Milky Way 


The number of known super star clusters in the Milky Way is 13, including RCW 38 


and the 12 clusters listed in Table 2 of 


Portegies Zwart et al. 


(120101) . The cluster members 


of RCW 38 were previously estim ated to be a few 1000 
was not included in the table by 


Portegies Zwart et al. 


Wolk et al. 2006) and RCW38 


(120101) which limited the cluster 


mass to be more 


stars flKnhn et al 


hanMO'^MQ. The recent estimate shows that RCW38 harbors 10,000 


20151) and it is appropriate to include RCW 38 as one of the super star 


cluster s. Among t 


1 Myr (I Wolk et al 


re 13 RCW 38 is the youngest one with an estimated age of less than 


2006 


, Sections 4.2 and 4.3 of the present work) and the ages of the others 


ar e > 2 Myrs. The Orion Nebul a Cluster is listed as an OB association with an age of 1 Myr 


m 


Portegies Zwart et al. 


(120101) . The ambient interstellar medium will reflect youth of a 


cluster, and cloud dispersal by the feedback will be less intensive in younger clusters than in 
more evolved clusters. According to eye inspection of the infrared images, only hve clusters 
among the 13 clusters, RCW38, NGC3603, Westerlund2, DBS[2003]179, and Trumplerl4, 
are associated with dust nebulosity, and their ages are from 2 Myrs to 3.5 Myrs except 
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for RCW 38. The remaining 8 clusters, most of which have ages between 3.5Myrs and 
18Myrs, have no associated nebulosity, with an exception of the Arches cluster having an 
age of 2 Myrs. It is likely that the interstellar medium around the 8 clusters is completely 
dispersed by the stellar feedback including ionization and stellar winds. 

Figure [H] shows the averaged CO intensity in the three clusters with nebulosity 
RCW 38, NGC3603, and Westerlund2, and indicates that RCW 38 still has signihcant 
molecular mass within ~3pc of the cluster. In Westerlund2 and NGC3603 whose ages 
are 2 Myrs the central part of the molecular clouds within 10 pc of the cluster is already 
strongly ionized/dispersed. The rich remaining gas toward the cluster is a unique feature 
of RCW 38 and is consistent with its age signihcantly smaller than 2 Myrs. RCW 38 may 
therefore provide a best opportunity to study the initial conditions for the super star cluster 
formation. 


4.2. The two clouds; physical association and evidence for collision 

The present observations have revealed the distribution of the two molecular clouds 
with I2kms“^ velocity separation toward RCW 38 at a spatial resolution of ~0.1pc. The 
following pieces of evidence offer robust verihcation of the association between the two 
clouds and the cluster. 

1. The ring cloud shows a central cavity created by the feedback of the cluster. The ring 
cloud is also associated with the infrared hlamentary/bubble-like features heated by 
the cluster. The huger cloud shows good correspondence with IRS I and the western 
cavity. The morphological correspondence supports the physical association of the 
clouds with the cluster. 

2. The line intensity ratio of the ^^CO J=3-2 emission to ^^CO J=I-0 emission shows 
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significantly high values, indicating high gas temperatures toward the cluster. LVG 
calculations show that the temperature there is higher than 20-30K, consistent with 
local heating by the O stars in the cluster. 

3. Another sign of the connection between the two clouds is seen as the bridging features 
in velocity between them at least in three places including the direction of the cluster 
center. Alternative possibility that the bridging features are due to stellar winds is 
unlikely as discussed later in this section. 


previous two cases Westerlund2 and NGC3603 flFurukawa et ah 


.2kms ^ as seen in the 

2009; 

Fukui et al. 

— ^ 
r—1 

o 

CM 


The velocity is too large for the clouds to be gravitationally bound by the total mass of the 
clouds and cluster, less than 10® M©. This indicates that the association of the two clouds 
is by chance. We frame a hypothesis that the two clouds collided with each other recently 
and that the collision triggered formation of stars, primarily O stars, as have been suggested 
for Westerlund2 and NGG3603. We note that the fourth and hfth super star clusters with 
nebulosity, DBS [2003] 179 and Trumplerld, are also associated with two molecular clouds 
which appear colliding to form the cluster (Kuwahara et ah in preparation; Fukui et al. 
in preparation). Therefore, it is possible that the hve super star clusters with nebulosity 
are all formed by triggering under cloud-cloud collision, suggesting the important role of 
collision in formation of a super star cluster and multiple O stars. In the remaining 8 super 
star clusters without associated nebulosity it is not likely that we are able to detect the 
parent molecular clouds. 


Stellar wind accel eration may be an alt ernative interpretation to cloud-cloud collision 


(e.g., see discussion in 


Furukawa et al. 


20091) . We see a velocity gradient in the huger cloud 


of RGW 38 (Section 3), which may suggest some acceleration by the O stars via the stellar 
winds. This is however not likely because the gradient is extended outside the cluster and 
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rather uniform over the hnger cloud; the velocity shift should be most prominent toward the 
cluster if the winds originating from the O stars play a role. We also note that the hnger 
cloud is located on the nearside of the cluster as suggested by the coincidence with IRS 1 
fFigure fT3|) . implying that the wind acceleration will cause blue-shift instead of the observed 
red-shift. The collision between the two clouds is therefore a plausible interpretation, and 
we exclude wind acceleration as the cause of the two velocities of the clouds. 


The possible role of the winds in acceleration raises an alternative to the origin of the 
bridging feature toward the cluster center. Since the bridging feature, having a tight link 
with the hnger cloud, shows a good correspondence with the O stars (Figure [TT]i), the 
bridging feature may be part of the hnger cloud accelerated by the stellar winds and not 
created by the collision. The mass of the bridging feature in Figure lid is estimated to be 
17 Mq from the intensity with the X(CO) factor. The kinetic energy in the bridging 

feature is estimated to be 10^® erg, small enough to be accele rated by the winds of ~20 O 


stars whose kinetic energy can be more than 10^^ erg (e.g., in 


Furukawa et ah 


2009 1. In this 


scenario, the wind-accelerated gas will be blue-shifted since the hnger cloud lies in front 
of the O stars. This contradicts red-shift of the bridging feature. The velocity range of 
the bridging feature which lies only between the two velocities of the collidi ng clouds, not 


beyon dt 


2015alJb : 


le velocity of t 


Takahira et al. 


le hnger cloud, is also consistent with the collision flHaworth et al. 


2014J: 


Torii et al 


201511 . The other two bridging features in the 


south (Figure [6]) are not associated with O stars and it is unlikely that the winds are a 
source of acceleration. To summarize, the three bridging features are likely formed by 
the collisional interaction between the ring cloud and the hnger cloud, lending additional 
support for the cloud-cloud collision. 


Recent theoretical studies of global numerical simulations of galactic disks indicate that 
collisions between giant molecular clouds are frequent with a mean free interval of ~10 Myrs 
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(Tasker & Tan 2009 

Dobbs et al. 2015 

). Th 

is is comparable to the tvpical lifetime of a giant 

molecular cloud 201' 

dvrs 

iFnkni et al. 

1999; 

Yamaguchi et al. 2001; Kawamura et al. 2009: 

Fukui & Kawamura 

2010' 

, suggesting that a giant molecular cloud experiences cloud-cloud 


collision more than once in its lifetime. It is also possible t hat the relatively rich m olecular 


gas in the Vela molecular ridge where RCW 38 is located flYamagnchi et al 


1999) favors 


more frequent encounters between the clouds than the average over the Wilky Way. 


4.3. Details of the collisional interaction 


Currently, we see the two clouds, the bridging feature, and the O star cluster, all of 
which are within 1 pc of the cluster center. The ring cloud has a cavity of 0.3 pc radius 
which is likely ionized by IRS 2 and the other O stars. The huger cloud is probably being 
ionized and observed as IRS 1 whose surface is irradiated by the O stars. Theoretical studies 
of the collisi on scenario show that a shock-compresse d interface layer is formed between the 


two clouds (iHabe fc Ohta 


1992 : 


Anathpindika 


20101) . It is likely that major part of the 


interface layer has already been converted to the O stars and that the bridging feature is a 
remnant of the interface layer. 

By assuming tentatively that the relative cloud motion is inclined to the line of sight 
by 45°, we estimate the relative velocity of the two clouds to be 17kms“^. The actual 


velocity of the shocked interface layer 
collision between the two clouds (e.g. 


may be decelerated 


Haworth et al. 


Dj momentum conversation in 
2015aJ), and we adopt 10 kms“^ as 


the relative velocity of the interface layer. The distribution of the triggering-formed O stars 
is extended by about 1 pc (Figure[II]). We estimate the typical collision timescale to be 
~ 1 pc/10kms“^ ~ 0.1 Myrs, which gives roughly the age of the O stars in the cluster. 


If we assume the O stars are formed by a constant mass accretion rate in the collision 
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time scale 10^ yrs, the mass accretion rate is estimated to be 4 x lO“^M 0 /yr for a 40 Mq 
star (an 05.5 star). Magnetohydrodynamical numerical simulations of cloud-cloud collision 
have shown th at the turbulence is e nhanced and the magnetic held is amplihed in the 


interface layer fjinoue fc Fukui 


20131 ). This leads to two orders of magnitude increase of 


the mass accretion rate from the pre-collision state; e.g., a core of 126 Mq is formed in 
0.3Myrs in the shocked interface layer as shown by model 4 of Inoue and Fukui (2013) and 
such a massive core is a promising precursor of an O star. The total mass of the O stars in 
RCW 38 is estimated to be ~ 20 x 20 = 400 Mq, if we assume a typical mass of an O star 
to be 20 Mq. This corresponds to ~10% of the total clump mass of the ring cloud 

3000 Mq (see Tablet). 

The cavity of the ring cloud is created by ionization of the O stars. The propagating 
velocity of the ionization front is roughly estimated to be 3 kms“^ from a ratio 0.3pc/10^yrs. 
This is consistent with the ioniza tion front velocity in a molecular cloud driven by Hll 


regions fjSpitzer fc Tomasko 


19681) . In Figureim RCW 38 shows a molecular peak within 


3 pc of the center, whereas the other two NGC3603 and Westerlund2 show decrease of 
molecular gas within 10 pc of the cluster. This trend is explained as due to ionization. The 
latter two clusters have an age of 2Myrs, an order of magnitude larger than RCW 38. By 
taking a ratio of 10 pc/2 Myrs, we roughly estimate the average velocity of the ionization 
front to be 5kms“^, which is comparable to 3kms“^ in RCW 38. If we assume velocity 
of the ionization front to be 3-5kms“^ and common molecular cloud density, the ring 
cloud of Ipc radius will be ionized within 0.2-0.3 Myrs and triggered star formation will 
be completely terminated. This poses an upper limit of a few times 0.1 Myrs for the age 
spread of the stars formed by triggering. We note that the age spreads of the youngest 
members of the two super star clusters, NGC3603 and Westerlundl, are determined to 
be 0 . 1 -0.3Myrs, an order of magnitude larger than their cluster age 2-5Myrs, by a 
careful photometric study with VLT for proper-motion selected young cluster members 
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with HST flKudrvaytseva et al. 


Baneriee fc Kroupal (120131) and 


2012). We note that moc 


el calculations of NGC3603 by 


Baneriee fc Kroupal (120141) also suggest the rapid cluster 


formation through a single starburst event followed by signihcant residual gas explusion. 
The small age spread is possibly explained by the collision scenario so that formation of ~20 
O stars by triggering is terminated quickly due to the ionization in the order of O.lMyrs. 
Most of the low-mass stars were probably formed in the ring cloud prior to the cloud-cloud 
collision, since they are extended over the ring cloud with a radius of a few pc (Figured]). 
We suggest that this naturally causes duality in the stellar age in RCW 38; i.e., young O 
stars of 0.1-Myr age and older low-mass stars of 1-Myr age, which is to be observationally 
conhrmed for the cluster members at infrared and X r ays. Such duality in ag e is also found 


in the other super star clusters (e.g., for NGC3603 see 


Haravama et ah 


20081 ). 


4.4. Conditions for the multiple O star formation 


It is shown that cloud-cloud collision triggered formation of the O stars in center of the 
super star cluster RGW 38. The present two clouds show a sign of collisional interaction, 
the bridging features, in the three regions denoted as BRl, BR2 and BR3 in Figure5. 
Only BRl shows formation of multiple O stars with projected density of 30(pc)-2 
[= 15/(1.Opc X 0.5pc)], and BR2 and BR3 no O star formation. The ring cloud has the 
highest molecular column density 1 x 10^^ cm“^ around BRl and has lower column density 
of (1 — 2) X lO^^cm”^ in BR2 and BR3. On the other hand, the huger cloud has nearly 
uniform column density of 1 x 10 ^^ cm“^, and the relative velocities between the two clouds 
are almost the same in the three regions. We suggest that the high column density is 
the necessary condition for O star formation in cloud-cloud collision. The total molecular 
masses toward the three regions are estimated to be ~260 Mq, ~ 9 OM 0 , and ~ 1 OOM 0 , 
respectively, suggesting that BR2 and BR3 are probably not massive enough to form 
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multiple O stars like in BRl, while it is possible that the collision is in its early phase prior 
to O star formation in BR2 and BR3. 


Tabled compares the observational properti es of the molecu 


the other two s uper star clusters, Westerlund2 (iFurukawa et ah 


ar clouds in RCW 38 with 


2009 


and NGC3603 


and RCW 120 flTorii et al 


Fukui et a. 


20 


201i 


Ohama et ah 


and with the single O stars in M20 flTorii et ah 


20101 


2m, 


), where cloud-cloud collision is suggested to be a trigger 


of O star formation. The comparison shows a trend that the multiple O star formation 
as in RCW 38 is found only for a high column density of 10^^ cm“^ in either of the two 
clouds, whereas the cloud mass and relative velocity do not show a particular trend in 
multiple O star formation; i.e., the mass range is from 10^ Mq to 10^ Mq and the observed 
relative velocity is 10-30 km s“^ in either of the multiple or single O star formation. The 
high column density in one of the two colliding clouds is therefore possibly the necessary 
condition for the multiple O star formation. 

The two scenarios for high-mass star formation, the monolithic collapse and 
competitive accr etion, are not yet verihed by confrontation with observations 


(ITan et al. 


20141) . Cloud-cloud collision, another possible scenario for high-mass 


star formation, provides directly testable three observational signatures including 

(1) two clouds associated with young O star(s) with (2) relative velocity greater than 10 kms“^, 

and (3) the bridging feature between the two clouds, and RCW 38 show all the signatures. 


According to 


Inone fc Fnknil (120131 ). the effective sound speed in the shocked interface 


layer is isotropic and given roughly as the velocity separation of the two colliding gas. 

It is then possible that low velocity collision of molecular gas may lead to a lower mass 
accretion rate and formation of a lower mass star. Observations of low-mass star formation 


without 0 stars triggered bv cloud-clouc 

collision are presentee 

bv several authors 

(Duarte-Cabral et al. 

2010; 

Higuchi et al. 

2010; 

Nakamura et al. 

2012 

). The relative 
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velocity of these collisions is as small as a few kms“^, and leaves room for non-triggered 
star formation under a gravitationally bound system.lt remains to be explored what is the 
role of cloud-cloud collision in low-mass star formation. 


5. Conclusions 

We have carried out CO J=l-0 and J=3-2 observations toward the super star cluster 
RCW 38 with Mopra, ASTE and NANTEN2 mm/sub-mm telescopes. The main conclusions 
of the present study are summarized as follows; 

1. We have revealed distributions of two molecular clouds at velocities of 2kms“^ and 

14kms“^ toward RCW 38. The ring cloud (the 2kms“^ cloud) shows a ring-like 
shape with a cavity ionized by the cluster and has a high molecular column density of 
~ 10^^cm“^. The other, the huger cloud (the 14kms“^ cloud), has a tip toward the 
cluster and has a lower molecular column density of ~ lO^^cm”^. The total masses of 
the ring cloud and the huger cloud are 3.0 x and 1.9 x 10^ Mq, respectively. 

2. It is likely that the two clouds are physically associated with the cluster as verihed by 
the high intensity ratio of the J=3-2 transition to the J=l-0 transition, R 3 _ 2 /i_o, 
toward the cluster. The observed ratio indicates kinetic temperature of > 20-30K 
according to the LVG calculations, which is signihcantly higher than the canonical 
molecular cloud temperature 10 K with no extra heat source. The heating is mainly 
due to the O stars. The association is further supported by the distribution of the 
molecular clouds corresponding to the cluster and the infrared dust features; the 
ring cloud shows the cavity toward the cluster center and good coincidence with the 
extended dust features, and the huger cloud corresponds to the infrared ridge IRSl. 
In addition, the two clouds are linked with each other by bridging features in velocity 
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at least in three places including the direction of the cluster, supporting further the 
physical connection between the two clouds. 

3. The total mass of the clouds and the cluster is less than 10® Mq. The velocity 
separation of 12kms“^ is too large for the clouds to be gravitationally bound, and we 
suggest that the clouds encountered by chance. We present an interpretation that the 
two clouds collided O.lMyrs ago at 17 kms“^ with each other to trigger formation of 
the O stars in the cluster. We argue that the double 05.5 star (IRS 2) was formed 
within this timescale at an average mass accretion rate ~ 4 x 10“‘^MQyr“^. 

4. The tip of the huger cloud and the bridging feature connecting the two clouds well 
coincide with the O stars within 0.5 pc of the cluster center, indicating that the 
triggering happened only toward the inner 0.5 pc and the other member low-mass stars 
outside 0.5 pc are pre-existent, being extended over the ring cloud where collisional 
interaction is not taking place. ROW 38 is the third super star cluster alongside of 
Westerlund2 and NGC3603 where cloud-cloud collision triggered O star formation 
in a super star cluster, lending further support for an important role of supersonic 
collision in O star formation. Among the three super star clusters, RCW 38 is unique 
because it is the youngest cluster where the initial conditions prior to the 0-star 
formation still hold without signihcant cloud dispersal by stellar feedbacks beyond 

1 pc. 

5. The present hndings suggest a possibility to prepare a recipe for O star formation in a 
super star cluster. Multiple O star formation is triggered where collision of two clouds 
at velocity 10-30 km s“^ takes place. One of the clouds has a high column density 
of ~ 10^®cm“^, while the other has a lower column density. Since O star formation 
takes place only at the spots of the collision, the distribution of O stars keeps the 
distribution of the collisional interaction as long as the gravitational rearrangement 
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of O stars after the collision does not play a role. This picture provides a natural 
set up for the highly turbulent interstella r medium as the in i tial co ndition for O star 
formation along the line of suggestion bv IZinnecker fc Yorkel (120071) . 


In order to better establish the role of cloud-cloud collision in O star formation, we 
need an extensive systematic study of molecular clouds toward a number of young Hll 
regions and young stars whose mass exceeds 20 Mq. Such a study is on-going and will allow 
us to test how cloud-cloud collision triggers and regulates O star formation. It is imperative 
to attain very high angular resolution in imaging the molecular gas to resolve dense cores 
in the shocked region and ALMA is the most important instrument toward the goal. It is 
also a challenge to test if mini-starbursts in galaxies are triggered by supersonic collisions 
between molecular clouds. 
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The 3.6 /im image of RCW 38 obtained with the Spi tzer/IRAC o bserva tions 


2008h. YSOs and O star candidates obtained with 


Wolk et ah 
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(1201 ll) are plotted in (b). Red circles and white circles indicate the class 


O/I and flat spectrum YSOs, respectively, while black dots indicate the class II and III 
YSOs. White crosses indicate the can didate O stms. fb ) A close-up of the central region of 


RCW 38 from the VLT observations flWolk et ah 


20061) . Z band is shown in blue, H band 


is green, and K band is red. The bright infrared emission IRS 1 and IRS 2 are shown by 


arrows. 
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Fig. 2.— J=l-0 integrated intensity distributions of the ring cloud (a) and the huger 

cloud (b). The area presented in Figure[3] with the Mopra data is shown by black lines in 
(b). (c) Latitude-velocity diagram of the J=l-0 emission. Dashed line indicates the 

position of RCW 38 IRS 2, while thick dashed lines show the regions used for the molecular 


mass estimates. 
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Fig. 3.— Integrated intensity distributions of the ring cloud and the huger cloud are shown 
in the left and central four panels, and the Dec.-velocity diagrams are in the right two panels. 
Upper three panels show the J=l-0 emission, and lower three shows the J=l-0 

emission. In (a), (b), (d), and (e), crosses show the position of IRS2, while the vertical 
dashed lines indicate the integration range in the declination-velocity diagram in (c) and (f). 
Thin dashed lines indicate the observed area with Mopra. In (c) and (f) the vertical dashed 
lines indicate the integration ranges in the integrated intensity maps in (a), (b), (d), and (e). 
Thick horizontal dashed lines indicate the declination of IRS 2. Three arrows in (c) indicate 
the declinations of the bridge features, BRl, BR2, and BR3, shown in Figure|6l 
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Fig. 


4.— R.A.-velocity diagram of J=l-0 is shown. Integration range in Dec. is from 
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Fig. 5.— Enlarged integrated intensity distributions of the ring cloud in J=l-0 (a) 

and J=l-0 for a velocity range of —8.3-+9.1 kms“^. Crosses indicate the position of 
IRS 2. 
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Table 1: Parameters of the clumps I-IV 



R.A. 

Dec. 

wlsr 

Av 

size 

Mlte 

Name 

I'h m Sj 

r'"] 

[kms“^] 

[kms“^] 

[pc] 

IMa] 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

clump I 

8:59:9.2 

-47:29:1.5 

5.0 

4.2 

0.5 

1.0 X 10^ 

clump II 

8:59:12.5 

-47:30:16.9 

1.7 

4.1 

0.4 

8.8 X 10^ 

clump III 

8:59:8.4 

-47:31:29.7 

2.3 

5.3 

0.3 

5.7 X 102 

clump IV 

8:59:1.2 

-47:30:27.3 

2.3 

4.8 

0.4 

5.6 X 102 


Note. —Column (1): Name. (2,3): Position in J2000. (4, 5): Peak ulsr and velocity width Au (FWHM) 
of the J=l-0 profile, where fitting with a Gaussian function is used to derive the parameters. (6): Size 
of the clump, measured as geometric average of major and minor diameters at the integrated 

intensity 45Kkms~^. (7) Molecular mass within the clump size. 
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Fig. 6.— A comparison of the J=l-0 distributions between the ring cloud (color) 

and the hnger cloud (contours) is presented. BRl, BR2, and BR3 show the positions of the 
bridging features by dashed circles. 
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Fig. 7.— The Dec.-velocity diagram of the J=3-2/J=l-0 ratio. The integration range 
in the R.A. is the same as Figure[3](c). The dashed line indicates the position of IRS 2. 
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Fig. 8.— Curves of the J=3-2/J=l-0 ratio as a function of for various n(H 2 ) 

ranges estimated with LVG for X(CO)/{dv/dr) of 10“^ (kms“^ 
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Fig. 9.— LVG results on the n(H 2 )-Tk plane for the regions A-E in the ring cloud are shown. 

J=3-2/J=l-0 is plotted in black, and ^^CO/^^CO J=l-0 is in red. The regions A-E 
are plotted with black crosses in the CO integrated intensity map in the upper-left panel, 
where J=3-2 is shown in color, while J=l-0 is shown with contours. White 

cross indicates the position of IRS2. Spectra of J=l-0 (black), J=3-2 (red), 

and J=l-0 (blue) are also shown. Intensities of the J=l-0 spectra are multiplied 
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10.— LVG results for the regions F-J in the huger cloud in the same manner as Figure 9. 
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Fig. 11.— Comparisons between the two molecular clouds (contours) and the Spitzer 3.6 pm 
image. The velocity ranges for the left two panels (a) and (c) and for the right two panels (b) 
and (d) are —8.5—1-9.1 kms“^ and -1-9.1—1-17.5kms“^, respectively. White contours in (d) 
indicate the bridge feature BRl integrated over a velocity range -1-9.1-11.5 km s“^. YSOs 
and O star candidates are plotted with the same manner as in Figure[U 
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Fig. 12.— Velocity channel maps of the J=3-2 emission (contours) superimposed on 

the Spitzer 3.6 pm image. In the left-top panel, approximate distributions of the infrared 
hlaments (hlaments a-h) are depicted by the dashed lines. The parenthetic numbers indicate 
peak velocities of the hlaments in kms“^. (see Table 2) 
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Table 2: CO velocity range corresponding to the infrared filaments in Figure[T2] 


Name 

(1) 

wlsr range 

[kms“^] 

(2) 

a 

-1.7~+0.7 

b 

-1.7-+0.7 

c 

+0.7-+3.1 

d (two filaments) 

-1.7-+0.7 

e 

+0.7-+3.1 

f 

-0.7-+3.1 

g 

-1.7-+0.7 

h 

1 

1 

1 


Note. — Column: (1) Name of the filament. (2) The velocity range in which CO emission corresponds to 


the filament. 
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Fig. 13.— J=3-2 distributions of the ring cloud (blue contours) and the huger cloud 

(red contours) superimposed on the VLT image in Figure[T] The contours are plotted at every 
15Kkms“^ from 130Kkms“^ for the ring cloud and at every 5Kkms“^ from 20Kkms“^ 
for the huger cloud. 
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Fig. 14.— Radial distributions of the CO emission in the three SSCs RCW38 (a), Wester- 
lund2 (b), and NGC3603 (c). (a) The NANTEN2 CO J=l-0 data set shown in Figure[2]is 
used. The present velocity range i n clude s the two colliding clouds, (b) The NANTEN2 CO 


J=2-l dataset in 


Furukawa et al. 


(120091) is used. Here the present velocity range includes 


only the blue-shifted cloud, because the red-shifted cloud which overlaps the SSC along the 
line-of-sight is indeed located rela tively distant from the cluster compared wi th the blue- 


shifted cloud (jOhama et al. 


2010). (c) The NANTEN2 CO J=2-l dataset in 


Fukui et al. 


(1201411 is used. The velocity range includes the two colliding clouds. 
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Table 3: Comparison of five regions of cloud-cloud collision. 


Name 

(1) 

cloud masses 

[Mq] 

(2) 

column densities 

[cm“^] 

(3) 

velocity separation 

[kms“3] 

(4) 

# of 0 stars 

(5) 

reference 

(6) 

RCW38 

(2 X 10^, 3 X 10^) 

(1 X 10^3, 1 X lO^b 

12 

- 20 

This study 

NGC 3603 

(7 X 10b 1 X 10b 

(1 X 10^3^ 1 X 1022) 

20 

- 30 

[1] 

Westerlund 2 

(8 X 10b 9 X 10b 

(2 X 10^3^ 2 X lO^b 

13 

14 

[2, 3] 

M20 

(1 X 10b 1 X 10b 

(1 X lO^b 1 X 10^2) 

7 

1 

[4] 

RCW 120 

(4 X 10b 5 X 10^) 

(8 X lO^b 3 X 1022) 

20 

1 

[5] 


Note. — Column: (1) Name. (2, 3) Molecular masses and column densities of the two colliding clouds, 
(blue-shifted cloud, red-shifted cloud) (4) Relative radial velocity between th e t wo clouds. (5) Number of O 


stars created via clo ud-c loud collision. (6) Ref erences: [1] 


Ohamaetaf 


(I 2 OI 0 I) . [4] 


Torii et al 


(I 2 OI 1 I), [5] 


Torii et al 


Fukui et al. 


(l2015h . 


(I2OI4I) . [ 2 ] 


Furukawa et al. 


(120091), [3] 






























